ASYMPTOTICS FOR PRODUCTS OF CHARACTERISTIC POLYNOMIALS IN 

CLASSICAL /3-ENSEMBLES 
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Abstract. We study the local properties of eigenvalues for the Hermite (Gaussian), Lagucrre (Chiral) 
and Jacobi /3-enscmbles oi N X N random matrices. More specifically, we calculate scaling limits of the 
expectation value of products of characteristic polynomials as N oo. In the bulk of the spectrum of 
each /3-ensemble, the same scaling limit is found to be e^^iFi whose exact expansion in terms of Jack 
polynomials is well known. The scaling limit at the soft edge of the spectrum for the Hermite and Laguerre 
/3-ensembles is shown to be a multivariate Airy function, which is defined as a generalized Kontsevich 
integral. As corollaries, when /3 is even, scaling limits of the fc-point correlation functions for the three 
ensembles are obtained. The asymptotics of the multivariate Airy function for large and small arguments 
is also given. All the asymptotic results rely on a generalization of Watson's lemma and the steepest 
descent method for integrals of Selberg type. 
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1. Introduction 

1.1. /3-Ensembles of random matrices. In this article wc consider three classical /3-cnscmbles of Ran- 
dom Matrix Theory, namely the Hermite (Gaussian), Laguerre (Chiral), and Jacobi /3-ensembles (II/3E, 
L/?E, and J/3E for short). Their eigenvalue probability density functions are equal to 

Xkf, X. eM, H/3E, (1.1) 
Xkf, e]R+, L/?E, (1.2) 

X, e(0,l), J/3E. (1.3) 

The normalization constants are all special cases of Selberg's celebrated formula |23| and are given in the 
appendix. 

For special values of the Dyson index j3, we recover more conventional random matrix ensembles |23| 
|4T| . The /? = 1, 2, 4-ensembles indeed correspond to the ensembles of random matrices whose respective 
probability measures exhibit orthogonal, unitary or symplectic symmetry. 

For general /3 > 0, Dumitriu and Edelman |13| constructed tri-diagonal real symmetric matrices with 
independent entries randomly drawn from some specific distributions and whose eigenvalues are distributed 
according to p.ip and (|1.2p . Killip and Ncnciu [3D] later obtained a similar construction for the J/3E. 
These explicit constructions play a key role in connecting the /3-ensembles with one-dimensional stochastic 
differential equations in the limit N ^ oo [T71I151I35] . Many probabilistic quantities of interest such as the 
global fluctuations, the gap probabilities, and the distribution of the largest eigenvalues were also studied 
in the limit N oo (sec for instance [8j[T5l[25,,29,46,.50j ) . More recently, some universality results (see 
below) concerning the general (3 > case have been obtained (see [19] and references therein) . It was shown 
for instance that as — oo, the eigenvalues in the bulk (middle) of the spectrum of any /3-ensemble are 
correlated, when appropriately rescaled, as the eigenvalues of the Hermite /3-ensemble. 

Apart from being related to the eigenvalues of random matrices, the densities (|l.ip ~ (|1.3p have alternative 
physical interpretations. Indeed, these densities appeared very recently in theoretical high energy physics 
[B, 9, 42, 48]. Moreover, densities such as (|l.ip ~ (|1.3p are equivalent to the Boltzmann factor for classical 
log-potential Coulomb gas and to the ground state wave functions squared for Calogero-Sutherland A'^- 
body quantum systems of the type A^v-i, and BCjy. We refer the reader to Forrester's lectures [35] 
for more details. The wave functions of the Calogero-Sutherland models are typically written in terms of 
a very special family of symmetric polynomials, namely the Jack polynomials [27[I37[|49] . This connection 
between the /3-enscmbles and the Jack polynomials has been exploited by many authors and has shown to 
be very fruitful [4,10 16,18,21,24,32,38 40,43] 

1.2. Products of characteristic polynomials. Now let X he an N x N random matrix in some (3- 
ensemble. Our aim is to find exact and explicit expressions for the large N limit of the expectation value 
of Y[j=i det(X — Sj). We will thus study the following expectation value: 

N n 

Kn{s,,...,s^)^{111[{x.~s,)) (1.4) 

1=1 j=i 

where x = {xi, . . . ,xn) denotes the eigenvalues of the random matrix X and where the angle brackets 
stand for the expectation value. More explicitly, for the densities ()l.l|) - ()1.3p . 

KNisi,- ■ ■ , s„) = 

N n p N 

X\X\{xt- Sj)cyi^^- -^V{xj)^ n \xi-Xjf dxi---dxN, (1.5) 

j=li=l 3 = 1 l<i<j<N 



G 



l<i<N 



-PXi/2 



n 



l<j<fe<Af 



n xf^e-^-/^ n 



l<j<fe<Af 



5a'(Ai,A2,/3/2) ^1 



<N 



i<j<k<N 
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where Zjsi is some normalization constant and 

f x] H/3E 
V[xj) = \xj~ (2//3)Ai Inxj L/3E (1.6) 

[-(2//3)Ailnx, - (2//3)A2ln(l-x,) L/3E 

Actually, we will see that it is more convenient to consider the weighted quantity 

(pjv(si,---,s„) =exp| - -^y(sj)| A'Ar(si,...,s„). (1.7) 

i=i 

At this point, it is worth stressing that if /3 is even and if we let n = A:/3, then ip^^si, . . . ,s„) gives 
access to the fc-point correlation function |23[BT| : 



Rk^Nixi, . . . ,Xfc) = 

k+N 

- I ... 

Nl Zk+N 



^ / • • • / cxp I - - ^ l^(a;_,)| ]^ |xi - Zjl'^dxfe+i • • -(ixfe+Ar (1.8 



j = l l<i<j<k+N 



Indeed, 



Rk,N{xi, . . . ,Xfc) = ^l^-*' Jl (X, -Xj)'^[v?Ar(si,...,S„)]^^j^^^j. (1.9) 

l<z<j</e 

Here the notation {s} > {x} means that the variables Si are evaluated as follows: 

S(i_i)/3+j = I'i for i = l,...,k and j = 1, . . . , ^. (1.10) 

For the special values /3 = 1, 2, 4, the scaling limits of products of characteristic polynomials are already 
well established (see [l][5j[7] and references therein). Thanks to the orthogonal polynomial method, they 
can be expressed as determinants or Pfaffians of one- variable special functions and their derivatives. The 
functions in question depend on the bulk or edge of the spectrum we are looking at |20| . Close to the hard 
edge of the spectrum (which correspond to Si ~ for Laguerre and s,; = 0, 1 for Jacobi ensembles), one 
gets a Bessel function or equivalently a o^^i (z) function. In the bulk of the spectrum of each ensemble (for 
instance, about Si = 1/2, 2A^, and for Jacobi, Laguerre, and Hermite, respectively), the formulas involve 
trigonometric functions or complex functions of exponential type, such as QFo{iz) ~ iFi{a]a;iz). Finally, 
at the soft edge of the spectrum (i.e., about sj = 4N for Laguerre and sj = ^/2N for Hermite) the scaling 
limits contain the Airy function Ai. In fact, these three regimes of large N asymptotics (Hard-Bcssel, 
Bulk- Trigonometric, Soft- Airy) correspond to the three most common universality classes for ensembles of 
random matrices with /? = 1,2,4 (e.g., see Chapter 7 in [23]). 

Much less is known about the general (3 > case. For the H/3E, Aomoto [3] and more recently Su j47) , 
obtained the limiting expectation value of the product of n = 2 characteristic polynomials respectively in 
the bulk and at the soft edge of the spectrum. When both n and N are finite but arbitrary. Baker and 
Forrester 2] proved that Kpf[si, . . . , s„) is either a multivariate Jacobi, Laguerre or Hermite polynomials 
with parameter a ~ [3/2 depending on whether the density considered is (|1.3p . (|1.2p or p.ip . They also 
used the theory of multivariate hypergeometric functions developed by Kaneko [55] and Yan [FT] , to express 
the expectation value Kj^(si, . . . , s„) in the Jacobi and Laguerre /3-ensembles as an n-dimensional integral: 

KN{si,...,Sn) = CN f ... f fle-^'P^'^^ n \tj-tk\^^^qN{t;s)dh---dtr.. (1.11) 

J'iS' •''^ j = l l<j<k<n 

where and qN{t] s) respectively denote the circle in the complex plane (or intervals of real numbers) and 
multivariate hypergeometric function whose precise forms depend on the ensemble. 

Such passage from an A'^-dimensional integral with Dyson parameter /? to an rt-dimensional integral with 
Dyson parameter /?' = 4//3 is an example of duality relation, which turns out to be very useful since the 
second n-dimensional integral allows us, in principle at least, to take the limit — > oo. As observed in [3], 
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when Si = • • • = s„, the function gAr(i; s) greatly simphfies. This was exploited in for determining the 
asymptotic behavior of the eigenvalue marginal density pm{x). 

Other dualities for the Hcrmitc and Lagucrrc /3-Enscmblcs were obtained in |10j . One particular duality 
was used to prove that as — > oo, for /3 — 1,2, 4, and for an appropriate choice of A and B, the expectation 
K]\j{A+Bsi, . . . , A+Bsn) in the H/3E is proportional to the folowing generalized Airy function (see Section 
[2] for more details about the notation): 

Ai('3/2)(si,...,s„) = -i- / e^P^^"'y^\Aiw)\^/%,4^^'\s,,...,s,,;tw)d"w, (1.12) 

which is absolutely convergent for all (si, . . . , s„) G M" and /3 e M+. The case /3 = 2 is proportional to 
Kontscvich's matrix Airy function |33] . 

1.3. Main results. We prove that at the soft edge, the expectation value of products of characteristic 
polynomials in both L/3E and H/3E actually lead to the same multivariate Airy function. Note that the 
simplest asymptotics for the multivariate Airy function are given in Proposition 13 . lOl 

Theorem 1.1 (Soft edge expectation). Let 

^(2^)V^ 2-V.A.-Va f^rtheWE, 

[AN, 2(27V)i/3 for the L/3E. ^ ^ 

Then as N ^ oo, 

^n{A + Bsi, ...,A + Bsn) - (27r)"(r^^,„)"' Ai^"/^\si, s„), (1-14) 
Note that the coefficients $jv,,i and r^^„ are given in the appendix. 

The above theorem suggests that the multivariate Airy function is the universal expectation value at the 
soft edge. In other words, for any /3-ensemble characterized by a potential V, the average of the product of 
n characteristic polynomials, when appropriately rescaled and re-centered at the soft edge, should become 
independent of V and should be proportional to Ai^'^/^-'(si, . . . , s„) as iV — ?> oo. 

The indication for universality is even stronger in the bulk of the spectrum. We indeed find that the 
three classical /3-ensembles possess the same asymptotic limit for the weighted expectation value ip^ in 
the bulk, which turns out to be a multivariate hypergeometric function of exponential type. We only state 
the result in the case of n = 2m, for simplicity; for n = 2m — 1, the combination oi ipN and (pN-i exhibits 
a universal pattern, which is given in Theorem 14. II 



Theorem 1.2 (Bulk expectation). For the H/3E, L/3E, and J/3E, let A be equal to \/2N, AN, and I, 
respectively. Let 



l^/Y^, ue{-\,\), H/3E, 
p{u)^{\4^ , ^e(0,l), L/3E, (1.15) 

Assume moreover that n = 2m is even. Then as N ^ oo, 

-^^^(Au+-^,...,Au+^) ^ 7™('V/?)e-^-P^(^) ,Fi^/'\2m/(3;2n/l3;2tns) (1.16) 

WiV,2m V f\u)N f\u)N / 

where ^iv,2m and jmWl^ respectively stand for the coefficient given in CO|) and (|XT3| . 

It's worth emphasizing that the universal coefficient ^„iWI^j when /? = 2, is conjectured to be closely 
related to the moments of the Riemann's C- function |5T| . 

The hard edge also involves a single hypergeometric series, which is o-P'i''^^^- The latter can be seen 
as a multivariate Bcssel function. We may thus surmise once more that the asymptotic expectation at the 
hard edge is universal. 
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Theorem 1.3 (Hard edge expectation). Let B he equal to N ^ and N ^ for the L/3E and J/3E, respectively. 
Then as N ^ oo, 

J—Kn{Bsu...,Bs„)^ o-Pi'^^'^((2//3)(Ai+n);-si,...,-s„). (1.17) 

t.N,n 

The coefficient ^N,n is given in (jA.lOp . 

Before explaining how wc prove these theorems, let us give some immediate consequences. As previously 
mentioned, if /3 is even and if ri = 2m = /3fc, then scaling limits of the fc-point correlation functions 
immediately follow from the above theorems. Since the hard-edge case for the L/3E and J/3E is already 
known (see Section 13.2.5 in [13]), we only display below the results for the soft edge and in the bulk. 

Corollary 1.4 (Soft edge correlations). Assume that (3 is even and n = (3k. Let A and B be as in Theorem 
\Lll Then as N ^ oo in the H/3E and L/3E, 

B'^Rk,NiA + Bxi,...,A + Bxk) ^ a^W) \A{x)f [A(^/'\s)] ^^^^^^^ , (1.18) 

where ak{(3) is given in (jA.lip . 

Corollary 1.5 (Bulk correlations). Assume that (3 is even and n ~ /3k. Let A and piu) be as in Theorem 
m Then as N ^ oo in the H/3E, L/3E, and J^E, 

A \k ( . Axi . Axk 



p{u)NJ "'"V p{u)N'"'' p{u)N 

bk{l3)\/^{2T:x)\f^ \e-'''P^'''\F[^''^\n/l3;2n/!3;2iT:s)\ , (1.19) 

where bk{j3) is given in (|A.12p . 

We stress that the function on the right-hand side of (|1.19p already appeared in Random Matrix Theory: 
it is exactly the same as the limiting fc-point correlation function of the circular /3-ensemble with (3 even, 
which is equal to the function p^/^^'^ixi, . . . , Xk) given in Proposition 13.2.3 of [53]. Moreover, very recently, 
the limiting fc-point correlation function in the bulk of the H/?E was shown to be universal (see Theorem 
1.2 in [H] ). Note that the latter reference does not give however, the explicit form of this universal 
fc-point correlation function. As a consequence of Corollarv ll.Si we now know that the universal fc-point 
correlation function, when (3 is even, is equal to the hypergeometric function on the right-hand side of 
(ICT . 

To the best of our knowledge, there is still no universality theorem regarding the limiting fc-point 
correlation function at the soft edge of the spectrum in /3-ensembles. However, assuming that such a 
universal correlation function exits for (3 even, we see from Corollary 11.41 that it must be equal to the 
fc- variable function on the right-hand side of (|1.18p . and as a consequence, it must involve the multivariate 
Airy function Ai*^^/^' . 

1.4. Organization of the article and proofs. As explained in Section[2l the multivariate hypergeomet- 
ric functions of the form pFq'^^ are defined as series of Jack symmetric polynomials. Although apparently 
complicated, these series can be evaluated efficiently by simple numerical methods |32] . 

The proof of Theorem 11.31 becomes trivial once we know that the expectations of products of charac- 
teristic polynomials in the L/3E and J/3E are respectively given by hypergeometric functions of the form 
^f['^''^^ and ^f['^''^\ This is known since [J. 

The proof of Theorems 11.11 and 11.21 is not so simple as that of Theorem 11.31 Indeed, the calculation of 
scaling limits in the bulk and at the soft edge requires the asymptotic evaluation of integrals of Selberg type, 
such as (|l.lip . The whole Section [3] is devoted to this task. We generalize the Laplace or steepest descent 
methods for higher dimensional integrals that contain the absolute value of a Vandermonde determinant. 

In Sectional we finally apply these results to our three classical /3-ensembles. Note that these explicit 
calculations abundantly make use of simple transformations of multivariate hypergeometric series, which 
will be given in the next section. Some remarks on PDEs satisfied by the scaling limits are given in Section 
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2. Jack polynomials and hypergeometric functions 

This section first provides a brief review of some aspects of symmetric polynomials and especially Jack 
polynomials. The classical references on the subject are Macdonald's book [37] and Stanley's article j49j . 
This will allow us to introduce the multivariate hypergeometric functions [28ll34l[5Tj . A few results proved 
here will be used later in the article. 



2.1. Partitions. A partition ti = K2, . . . , k^, . . .) is a sequence of non-negative integers Ki such that 

> K2 > • • • > > • • • 

and only a finite number of the terms are non-zero. The number of non-zero terms is referred to as 
the length of k, and is denoted 1{k). We shall not distinguish between two partitions that differ only by a 
string of zeros. The weight of a partition k is the sum 

|k| := ki + K2 + • • • 

of its parts, and its diagram is the set of points («, j) G such that 1 < j < k^. Reflection in the diagonal 
produces the conjugate partition k' = (k'j^, Kj, . . .). 

The set of all partitions of a given weight are partially ordered by the dominance order; k < a li and 
only if X^iLi '^i — SiLi '^i for all k. One easily verifies that n < a ii and only if a' < k' . 



2.2. Jack polynomials. Let A„(x) denote the algebra of symmetric polynomials in n variables a;i, . . . , x„ 
and with coefficients in the field F. In this article, F is assumed to be the field of rational functions in the 
parameter a. As a ring, A„(a;) is generated by the power-sums: 

Pk{x):=x\ + --- + xl. (2.1) 

The ring of symmetric polynomials is naturally graded: A„(a;) = ©fc>oAj;(a;), where A^j(a;) denotes the 
set of homogeneous polynomials of degree k. As a vector space, A^(a;) is equal to the span over F of all 
symmetric monomials 771^(2;), where k is a partition of weight k and 

mK{x) := Xi^ ■ ■ ■ x'^i" + distinct permutations. 

Note that if the length of the partition k, is larger than n, we set m^{x) ~ 0. 

The whole ring A„(x) is invariant under the action of homogeneous differential operators related to the 
Calogero-Sutherland models U : 

1=1 1=1 ' i<«#j<" 



The operators Ei and D2 are special since they also preserve each Aj"j(a;). They can be used to define the 
Jack polynomials. Indeed, for each partition k, there exists a unique symmetric polynomial pjf\x) that 
satisfies the two following conditions [49] : 

(1) Fi")(a;) =m«,(a;)-t-^c«^77v(a;) (triangularity) (2.3) 

(2) (^L>2 - ^(n - l)Ei^P^°''>{x) = e«P^")(2;) (eigenfunction) (2.4) 

where the coefficients and Ck^i belong to F. Because of the tiangularity condition, A„(a;) is also equal 
to the span over F of all Jack polynomials Pk"\x), with k a partition of length less or equal to n. 
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2.3. Hypergeometric series. Recall that the arm-lengths and leg-lengths of the box in the partition 
K are respectively given by 

anihj) = Hi" i and LH,]) = n'j - i- (2.5) 
We define the hook- length of a partition k as the following product: 

^«"^= n {i + a4hj) + ^Ui,j)). (2.6) 
Closely related is the following a-deformation of the Pochhammer symbol: 



n (--V) - n (-+<(^jvk(^,j)) (2.7) 



i<i<i(K.) ' (ij)eK 



In the middle of the last equation, {x)j = x{x+l) ■ ■ ■ (x+j — l) stands for the ordinary Pochhammer symbol, 

to which [x]!;"^ clearly reduces for i^n) = 1. The right-hand side of (|2.7|) involves the co-arm- lengths and 
co-leg-lengths box in the partition k, which are respectively defined as 

aUhi)^j-h and (2.8) 
We are now ready to give the precise definition of the hypergeometric scries used in the article. 

Definition 2.1. Fix p,q G = {0, 1,2,.. .} and let oi, . . . , Op, bi, . . . ,bq be complex numbers such that 
{i — l)/a — bj ^ No for alH € Nq. We then define the (p, g)-type hypergeometric series as follows [2811341151) : 



^FtHau....a,M,...,b,;x)=Y^ ^ ^^^T^^^^^P^'^'i-)- (2-9) 



Similarly the hypergeometric scries in two sets of n variables, x = {xi, . . . ,x„) and y ~ (j/i, . . . ,yn), is 
given in |51j by 



J^(")(a, a-b, b-x-y) - Y ' ■ ■ • [^.^ (y) 10) 



where we have used the shorthand notation 1" for 1, . . . , 1. 

Note that when p < q, the above series converge absolutely for all x € C", y € C" and a e R+. In the 
case where p = <?+ 1, then the series converge absolutely for all ||x|| < 1, \\y\\ < 1 and a e K.-|_. See [2H] for 
more details about convergence issues. 

Now we give some translation properties of o-^o"'' and iJ^q°'\ which prove to be of practical importance. 
For convenience, we write 

n 

^ ^ ^ / X\ Xjl \ 

(a") = (a, ...,a), b + ax ^ {b + axi, . . . ,b + aXn), = — , (2-11) 

1 — ax \ 1 — axi 1 — aXn ' 

where a, b are complex numbers and x = (xi, . . . , .t„). 
Proposition 2.2. We have 

oT^°'\a + x;b + y) ^ exp{nab + api{y) + bpi{x)} oT^"\x; y) (2-12) 

and 

n 

1 J-<"'(a; b + x;y) ^Uil^ by,)-^ i J"^"' (a; x; ^^). (2.13) 
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Proof. First set F{a, b) = QTlf\a + x; b + y) . Then, according to Eq. (3.3) of [4], 

E^>^\4"\x;y)=p,{x)„Tt\x;y), (2.14) 

so that 

\F. 



Similarly, 

dF 

— ^Pi{b + y)F. 
oa 

By solving those differential equations we get 

F{a,b) = e°fi(^+^'F(0,6) = e"fi(''+^'e^Pi(^)F(0, 0), 

which is the first desired result. 

For the second result, it suffices to prove 

n 

iJ-^"'(a;x;2/) - H)" = (a; x - 6; ^^). (2.15) 

Now let Gi{b) and Gr{b) respectively denote the left-hand and right-hand sides of (|2.15l) . 

In Eq. (A.l) of [1], we substitute x and y by cx and y/b, respectively. Then we let b,c ^ oo, and 
conclude that i!F^\a\x\y) satisfies 

E^^'^F - E^y^F = api{y)F. (2.16) 

From (|2.16l) . we get 

G'M - (-4'\^^"^(a;x;y)-}-i?<^\-F^"'(a;2;;y))|^^^_,^^^^ =-api(^)a(6). 

Finally, it's easy to check that Gi{b) also satisfies the differential equation just given above with the same 
initial condition at 6 = 0, which means Gi{b) = Gr{b). □ 

Corollary 2.3. 

oJ't\xi, • • ■ : a:n; a^ 6"-^-) = e'P'^-\F^''\k/a; n/a; (a - 6)xi, . . . , (a - b)xn) (2.17) 
= e('^-'')^-i+''fi(-)iF|")(fc/a; n/a; (a ~ b){x2 ~ x,), . . . , (a - 6)(x„ - x^)). (2.18) 
Proof. Proposition 12.21 implies that 

oJ-^^(xi, . . . , x„; a\ b"-") = e''Pi(-)oJ-<"'(a;i, . . . , a;„; (a - &)^ 0"-'=). 
The right-hand side is equal to 

^bM.) (a-6)l"lpi°^(x)Pi"\l^) 

i{K)<k V-"- J 

But we also know from Eq. (10.20) of Chapter VI [37] and Eq. (^77)) above, that 
Moreover, one easily checks that [fc/a]!"'' = whenever £{k) > k. Consequently, 

which is equivalent to the first equality. 

Likewise, the second equality also follows from Proposition 12.21 □ 
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We conclude this section with introduction of the fohowing symbol, for simplicity, 

Ef{x) Ef{x^, . . . , x„) = oJ-^"^(x; (-l)^ 1"-^), j = 0, 1, . . . , n. (2.19) 
These functions possess some similar properties with the exponential function, for instance, 

Ef{2i^ + x)^E'f\x). 

3. Asymptotic methods for integrals of Selberg type 

We want to get the large N asymptotic behaviour for integrals of the form (|l.lip by generalizing the 
classical steepest descent method, or more generally Laplace's method for contour integrals. We refer the 
reader to Olver's textbook [44] for more details on Laplace's method in the one-dimensional case. 

Before going any further, let us adopt some new notational conventions. Our variables are t = (ii , . . . , t„) 
and the Vandcrmonde determinant is given by 

n (3.1) 

l<i<j<Ti 

Moreover. A = (Ai,...,A„) denotes a sequence of parameters such that A^ > for all j while k = 
(ki, . . . , kn) denotes a sequence of non-negative integers of weight = kj. Finally, t^^^ = Y\ - 1^^ ^ 
and t'' = tf ■ 

3.1. Watson's Lemma. We first give the Watson's lemma for multiple integrals with Vandcrmonde 
determinants. The process of the proof will suggest a natural extension to contour integrals. 

Lemma 3.1. Let \\t\\ = max{|ti|, . . . , |t„|} and q{t) be a function of the positive real variables tj, such 
that 

m— 1 

'?W = *'"'(E«^-(^) + 0(11^11")) (II^II^O) (3-2) 

where 

|fc|=i 

is a homogenous polynomial of degree j . Then 

/•oo />oo m — 1 . 

In:=J^ ■■■j^ e-''^U^^\mfq[t)dt,---dt,, ^ N-J\x\+, +0[N--^-\>^\-^) (3.3) 

as N oo provided that this integral converges throughout its range for all sufficiently large N , where 
np ~ l3n{n — l)/2 and 

POO roo 

Aj^ •••/ t^"ie-^5'=i*^|A(t)|''aj-(t)dti---di„. (3.4) 

Proof. For each m, define 

m— 1 

0m(O 51 "jW- (3.5) 



One obtains 



I^-Y. + ■■■I e-''^U'^A{t)fc^,^{t)dt,---dt„. (3.6) 

As ||i|| — > we have (f>m{t) = i'^^"'^0(||t||™). This means that there exist positive constants k„i and K,n 
such that 

|</'™(t)| < A„||t||'"+I^l-" (0<||t||<fc„). 
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Accordingly, 



where 



e-^i:".i*.|A(t)|V™(i) dh---dt„ 



< 



e-i:;=it,|A(t)|'3ppl+™-"dti...d<„. 



(3.7) 



For the contribution from the range 

(O,oo)"\(0,fc,„)" =|J(Ji X ••• X J„), 

where Jj = (0,km) or [km, oo) but at least one of which is the infinite interval [km,oo). One knows there 
are (2" — 1) intervals in the union above. Without loss of generality, we only consider the case when 

J = (0, fc™)"-l X [k,n,00). 

Let Nq be a value of N for which the integral on the right-hand side of p.6p exists, and write 



so that $m(0 is continuous and bounded in J = (0,fc,„)" ^ x [fc,„,oo). Let i,„ denote the suprenium of 
|$,„(i)| in this range. When N > Nq, one finds by partial integration 



Ij:^ I er^^Ut^\A{t)\Um{t)dh---dtn 

g-(Af-JVo)Er=i*<$„^(i)lrfi,.^ . ..dU^dt„, (3.8 



j=0 l<»i<...<j,<„-l"'(0>'=^™)^x['^'"^°°) 

where [/(t)] denotes the evaluation of f{t) at ti = k„i if Z 7^ ii, . . . , 



Thus 



(0, /c x[A;^, 00) 



(3.9) 



Combining p.7p and p.9p . we immediately sec that the integral on the right-hand side of p.6p is 

0(7V-"pi-l->'l-™) as — > 00, and the lemma is proved. □ 

Corollary 3.2. With the same assumptions as in the previous theorem and with Xi — ^ > for all i, we 

have as N ^ 00, 



'(0,00) 
where 



e~JVE?=it.|A(t)|^g(i)d"i = Yl 

j=o 



r(l + /J/2 +j/3/2)r(/. + J/3/2) ao 



r(l + /3/2) 



0{N- 



-riyg —fin— 1 N 



«o = [t-^+'9(i)] 



t=(0,...,0) 



Remark 3.3. The integrals Aj, given by p.4p . is in general very difficult to evaluate. However, if q{t) is 
a symmetric function and Ai = ••• = A„, then Aj can be calculated by the Macdonald-KadcU-Selberg 
integrals, see [37|i P- 386 (a) or the original papers pTll^ . 
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3.2. Laplace's method: a single saddle point. Consider the integral 

In= •••/ cxp{-NY,Pitj)}h{A{t))q{t)dh---dt„, (3.10) 

J a J a j ^ 

in which the path ^ between the points a and b, say, is a contour in the complex plane C, p{x) and q(t) are 
analytic functions of x and t ~ {ti, . . . , tn) in the domains T C C and T", respectively. Here iV is a positive 
parameter and h{x) is a homogeneous analytic function of degree > 0, i.e., h{cx) = c'^h{x),\l Rec > 0. 

Like the one-dimensional integral, to obtain the asymptotics of integrals one usually needs to deform 
the path through some special points at which p' {x) = 0, called saddle points, we refer to sections 7 and 
10, |44j for more details about saddle point and paths of steepest descent. 

Recall that xq is a saddle point of order /x — 1 if 

p'(xo) = • • • = p(^-i'(.To) - 0, p(^)(a;o) ^ 0, 

where the integer /i > 2. In particular, when /i = 2 it is called a simple saddle point. The most 
common cases for the integrals are that p{x) has (1) one simple saddle point; (2) two simple saddle points; 
(3) one saddle point of order 2, at interior points of the integration path. Those occur ubiquitously in 
Random Matrix Theory, corresponding to the hard-edge, bulk and soft-edge limiting behavior. Fortunately, 
deformations to the path through saddle points for our multi-dimensional integrals can reduce to the 
one-dimcnsional case. But great care must be taken in dealing with the part involving Vandermonde 
determinant. 

We will first consider the case of a single saddle point of order — 1, which generalizes both cases (1) 
and (3). By convention ph(z) denotes the phase or argument of complex variable z. For an interior point 
xq of {a,h)gs, we denote 

w = angle of slope of ,^at xq = lim{ph(x — .tq)} (.t — x^ along (xq, h)^). (3-11) 

We moreover assume the following: 

(i) p{x) and q{t) are single-valued and holomorphic in the domains T C C and T"; h{x) is a homo- 
geneous analytic function of degree i/ > in C. 

(ii) The integration path is independent of N . The endpoints a and 6 of ^ arc finite or infinite, 
and {a,h)ss lies within T. 

(iii) p' [x) has exactly one zero of order /i — 1 at an interior point .tq of , i.e., p' [xq) = ■ . . = 
p(^"^)(a;o) = Q,p^'^\xo) ^ 0; in the neighborhoods of xq and t^ ~ (xp, . . . ,a;o), p{x) and q{t) can 
be expanded in convergent series of the form 

oo oo 

p(a;) =p(a;o)+^p.(a;-a:o)-^+'',g(i) = (i - io)^"' '^^ ~ 

where po ^ and q.j{t) = Y.\k\=3 Ij.k ■ 

(iv) There exists A^o > such that Ino converges absolutely at (a, . . . , a) and (&,...,&). 

(v) Re{p(x) — p{xq)} is positive on {a,b)^, except at zq, and is bounded away from zero as a; — > 
a or b along^. 

Remark 3.4. Let po := p^^^^ (xq) / l^l ■ When the phase of {x — xo)'^po is exactly equal to zero, then the first 
non-null term in Iie{p{x) — p{xq)} reaches its minimum value while Im{p(x) — p(xo)} gets equal to zero. 
Any path £P for x that guarantees ph{(a; — .To)^po} = is called the steepest descents. For our special cases 
coming from Random Matrix Theory, we will always be able to deform a part of the contour of integration 
to one of steepest descent contours. For the general case, however, the condition Kc{p{x) — p{xo)} > is 
sufficient. 

We proceed in three steps. 

Step 1. We introduce the following convention: the value of wq — ph(po) is not necessary the principal 
one, but is chosen to satisfy 

\UJ0+ ^iUj\<^TT, (3.12) 
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and the branch of ph(po) is used in constructing aU fractional powers of pq which occur. For example, Pq^^ 
means exp{(ln|po| + i'^o)//-*}- Since 

p{x) ~p{xo) ^ poix - xq)^, asx ^ xo along (a, b)^, 

introduce new variables wj by the equation 

= {p{tj)~PM)/Po- (3.13) 

The branch of ph(w_, ) is determined by 

/xph(wj) /io; (wj xq along (xo,b)^), (3-14) 

and by continuity elsewhere. 

For small ||t — ioll, Condition (iii) and the Binomial theorem yield 

Wj = {tj - Xo){l + -^{tj ~ Xo) ^ }. 

Application of the inversion theorem for analytic functions shows that for all sufficiently small p > 0, the 
disk \tj — xq\ < p is mapped conformally on a domain W containing Wj = 0. Moreover, if Wj G W then 
tj — xq can be expanded in a convergent series 

oo 

tj - Xo= y^csw-, 

s=l 

in which the coefficients Cs are expressible in terms of the ps- For example, ci = 1 and C2 = —pi/{ppo)- 

Let Ti,T2 be points of {a,xo)^, {xo,b)g/i respectively chosen sufficiently close to xq to ensure that the 
disk 

\w,\ < min{b(Ti) -p(xo)^/^ \p{r2)-p{xo)\'/''} 
is contained in W. Then {ti^xq)^ and {xq,T2)^ may be deformed to make its Wj map straight lines Ln 
and L211 respectively. If 

1^] = {p{tj) - P{xo))/po, j = 1, 2, 
then ill and L21 are directed line segments, respectively, from ki to and from to K2- On the other 
hand, let Li = L12 U Ln and L2 = ^21 U L22 denote the half-lines containing points 00, ki, and 0, K2, 00, 
respectively. Thus we have 

/ ••• / exp{-A^^p(t,)}/i(A(t))q(t)dti---dt„ = 

/ exp{-7Vpof^u;j^}/i(AH)/(u;)du;i---du;„, (3.15) 

J Lii\jL2i J Lii\jL2i 



-nNp(xQ) 



where 



/(-) = .(t) ( n E ^) ^ n ( E -.^-r^ 

i<j s>l * ^ j = l s>l 

For small /(w) has a convergent expansion of the form 



=ti;^-i^a,(«;), aj(?i;) - ^ aJ,fe^z;^ (3.16) 
]=a \k\=j 

in which the coefficients aj^k can be computed in terms of ps and qj^k- In particular, ao = go- 
Step 2. Following the approach for the Watson's lemma, we define fm{w), m = 0, 1, . . . , by 

m — 1 

f{w) = w^-^ J2 + w^^'fM. (3.17) 

3=0 
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Then fmiw) = 0(||w||"). Set = iyn{n - l)/2 and 

Aj — / ••• / exp{— po w^}h{A{w)) w^~^aj(w) dwi ■ ■ ■ dwn, (3.18) 

JLiUL2 JLiUL2 j—l 

the integral on the right-hand side of p.lSp is rearranged in the form 



/ •••/ cx.p{~Npo^w'^}h{A{w)) f{w)dwi---dw„ 

t/LiiUZ/21 '/LiiUZ/21 j—1 



A ■ 

E ^r(n.+m+.)/M - + (3.19) 



Here 

em,i{N) = E E / <^^v{~NpQY,w'^}KHw))w^-^aj{w)dwi---dwn, (3.20) 
summed over J = Ji x ■ ■ ■ x J„. = Lj^x U L21 or L12 U L221 but at least one of which is L12 U L22, and 



I ■■■ I cxp{~NpoJ2wj}HHw))w^-^UHdwi---dwn. (3.21) 

For Sm,2{N), splitting the integral domain into 2" parts, for one of which, for instance, the domain 
Lii X • • • X Lii X L21, substitute wj = KiVj,j < n and Wn = K2Vn, and note that Condition (v) implies 

Re{po<} > Vi, Re{po4} > Vi (3-22) 

for some positive r]i . we have 

p p ^ I*! I*! 

/ ••• / ex■p{~Npo^w'^}h{A{w))w^'^f„^{w)dwl■■■dwn^ ■■■ 
Jill JL21 Jo Jo 

exp{-Afpo(« + • • • + «_i + 4<}/i(A(«;)) ii;-^-iO(||u;|r") dvi ■ ■ ■ dv^ = 0(iV-("-+l^l+™)/'^). 

For em,i(-/V), without loss of generality, consider the ease when Ji = ■ ■ ■ = Jn-i = Lii U L21, Jn = 
Li2UL22- Then 

P n 

Ij : = / Gxp{—NpQ w^}h{A{w)) w^~^aj{w) dwi • • • dwn 

/n n 

exp{-(A^ - Nq)po E exp{-iVoPo E ^ ^^C^l^^)) w^^^aj{w) dwi ■ ■ ■ dwn- (3.23) 
j=i j=i 

Condition (v) implying that Re{poU'j } > 0, 1 < j < n for G Ln U L21, combining with p.22p and 
Condition (iv) we obtain 

|//| < C)(l)e-(^-^°^''i. (3.24) 
The combination of the results of this step with p.lSp leads to 

rT2 t'T2 ^ 

/ ••• / cxp{-Nj2pit3)}KMt))q{t)dh---dU = 

51 7v("-+I^I+.J')/a' +'^(Ar(n.+|A|+™)/A'''}' (^'^^^ 



m— 1 . 
-n7Vp(a;o) f \ " 



e . 

0=0 



as N ^ 00. 
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Step 3. It remains to consider the tail of the integral, that is, the contribution from (a, 6)" \ (n, T2)" = 
Uj Ji X • • • X J„, Jj ~ (ti, T2) or (a, ri) U (t2, 6) but at least one of which is (a, n) U (t2, 6). Without loss 
of generality, we assume that Ji = • • • = J„_i = (ti, T2) and J„ = (a, ti) U (t2, 6). From Condition (v) we 
know that 

Re{p(a;) -p(xo)} > 0(a; e (a, 6)^) and Re{p(a;) - p(xo)} > ?y2 (a; G (a, n)^ U (t2, &) (3.26) 
for sonic 772 > 0. Therefore, 

NKe{p{tn) - p{xq)} >{N- No)m + iVoRe{p(i„) - ^(2^0)}, 
and condition (iv) shows 



„-{N-No}i}2\„-nNop(xo) I 



/ exp{-iV ^(p(tj) - p(xo))}/i(A(i)) q{t) dti--- dU 

^ J A 1 

n 

/ \cxp{~NoJ2pitj)}HMt))<l{t)dti---dtJ. (3.27) 



i=i 

Thus the asymptotic expansion p.25p is unaffected by the tail integrals. 

By following the steps 1-3, we have established the following fundamental result 



Theorem 3.5. With the foregoing assumptions (i)-(v), Vm G N, we have 
j^ - J^ exp{-iV^p(t,)}MA(t))g(Odti---di„ = e-"^P(-»){^^^^^ 



- 0(- ' 



J=<J 

as N ^ 00. Here the coefficients Aj are given by (|3.18p . 

Remark 3.6. If we focus on the leading term with the coefficient Ag, then there is an immediate general- 
ization of the previous theorem which will be used later in the article. Let g{t) be analytic in the whole 
C". Now in the integral p.lOp . replace q{t) by q{t)g{N^^'^{t — to)), and assume /i > ^. At first sight, the 
factor g{N^^P'{t — to)) may become very large as — >■ 00. However, by repeating the analysis done in steps 
2 and 3 (and paying a particular attention to the rescahng of the variables), we see that giN^^f'it-to)) 
does not lead to a new significant contribution to the leading coefficient. Most importantly, if = /i, then 
we get a very simple formula: 



pb pb n 

/ •••/ exp{-NY,P{tj)}h{Mt))q{t)9{N^^''{t-to))dti---dtn 

•la J a 



-nNp(xo) f A'q ,r)(jyj-{n,, + \M + i)/f^)\ 



A'Q = qn ■ I exp{-po^w^}^(A(«;))it;^ ^ g{w) dwi ■ ■ ■ dw„ 

JLiUL-i JL1UL2 



where 

" ' ' (3.28) 

3.3. Laplace's method: two simple saddle points. Once again we consider the integral p.lOp . This 
time howerver, we suppose that p{x) has two simple saddle points x±, which means 

p'{x±) = and p± := p"{x±) ^ 0. 

This case is quite different from the one-dimensional case because of the Vandermondc determinant. Our 
assumptions are: 

(i) p{x) and q{t) are single-valued and holomorphic in the domains T C C and T"; h{x) is a homoge- 
neous analytic function of degree > in C; for simplicity, we also suppose that q{t) is symmetric 
and h{x) = h{—x). 

(ii) The integration path 0^ is independent of TV. The endpoints a and 6 of ^ are finite or infinite, 
and {a,h)<^ lies within T. 
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(iii) p'{x) has exactly two simple zeros at interior points x-^-, x- of setting tj^+ = (a;+, . . . , a;+, a;_ , . . . , x-) 
consisting of j x+'s and (n — j) x-'s, p± = p"{x±) and qj,+ = q{tj^+). 

(iv) There exists A^o > such that 1^^ converges absolutely at (a, . . . , a) and (&,..., 6). 

(v) Re{p(x) — p(a;ip)} is positive on {a,xo]a^ and [a;o,6),ga respectively for some xq G {a,b)^, except 
at x^^ x-i-, and is bounded away from zero as x — > a, or x — >■ 6 along . 

As in step 1 of Section [3?2l let Li,_|-Ui2,+ and Li^_UL2,- denote the unions of two half-lines respectively 
corresponding to 2;+ and x-. Thus, the n-dimensional integral can now be broken into a sum of 2" terms, 
each being an n-dimensional integral in which each variable is in the neighborhood of either x+ or x- . 
When j variables are in the neighborhood of x+ , the Vandermonde determinant becomes 

A,,+ H = (2;+ -x-p'--'') Yl {wp-w<j) Yl {wp-w<j). 

l<p<q<j j<P<9<" 

By the symmetry assumption, there are (") terms amongst the 2" terms that produce such a Vandermonde 
determinant. By proceeding as for Theorem [33] and focusing only on the dominant contribution we rapidly 
establish the following theorem which contains integrals similar to (|3.18p : 



/ expj - ^^wf - — ^ u;f |;i(Aj,+(w))dwi • --dw,,. (3.29) 



/(Li, + UL2, + )Jx(Li,_uL2,-)"-J' ;=i l=] + l 

Theorem 3.7. With the foregoing assumptions (i)-(v), 



/ ••• / cxp{~Nj2p{tj)}HMt))Qit)dti---dtn = 

•la J a j_2. 

(") cM-N{lp{x+) + (n - l)p{x^))} (qi^+ 5^^^7v-(''^+("-')-+")/2(l + 0(iV-i/2))) (3.30) 



E 

as iV —> 00. Here the coefficients Bi^^ are given by (|3.29p . 
Notice the fact that 

l„ + {n- = v{l - + vn{n - 2)/4 (? = 0, 1, . . . , n) 

attains its minimum value at I = m if n = 2m or I = m, m — 1 ii n = 2m — 1, which shows that probably 
only one or two terms in the sum of p.30p give a major contribution. 

3.4. Integrals of Selberg type. Recall that our aim is to evaluate the integrals such as when 
— > 00. In Theorems 13.51 and 13.71 h(A(t)) must be homogeneous and analytic. The latter condition 
is problematic since we want to calculate integrals involving |A(t)|'^. Of course, when t £ R" and v is 
even, then we can set h{A{t)) = (A(t))''. For ly not even, it will be enough to use the following integral 
representation: 

/>oo 

r-''-^H„{rx)dr, (3.31) 
which holds if x S R, > -1 and 7^ 0, 2, 4, . . .. Here 

^-W = E c. = -sin(^)r(i.+ l). (3.32) 

This can be computed by complex analysis method when — 1 < < 0, see Gradshteyn and Ryzhik's 
book |26| : otherwise, first integrate by parts and use the former result. 
Let us now consider 

In^ I ■ ■ I exp{-iVVp(tj)}|A(t)r<7(t)dti---di„, (3.33) 
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where the functions p and q are as mentioned previously while denotes some interval (a, 6) C M or one 
circle in the complex plane. Assume that v is not even and that the path = (a, b) C R, then 

In=c, f cxp{-NS2p{t,)}q{t) ( f r-''-^H,{rA{t))dr] d"t. 

J{a,b)^ \Jq ) 

By Fubini's theorem, we rewrite 



r r-^'-^i \ Gxp{-NY^v{t,)}q{t)H,{rA{t))d^t\dr. (3.34) 

"'0 Vj(a,6)" / 



In = c. 



Likewise, if v is not even and ^ is the unity circle T := {z : \z\ = 1}, setting tj = e'^^, since 
|A(t)| = | n (2 sin ^^) I and [] (2 sin ^^) = A(t) ^(^^,)-^"-l)/^ 

l<j<k<n l<j<k<n j=l 

we have 

./ r-^-^i exp{-Ny2p{t,)}q{t)H,{rA{t)T\{it,)-(^-^y^)d"t)dr. (3.35) 

Jo ViT" ^.^1 ^-^1 y 

Notice that the precise form of the integral over r in p.34p and p.35p allows a rescaling of the variable, 
so the function Hi,{r) plays a similar role like a homogeneous function of degree i^, and we can apply 
established theorems for the interior n-dimcnsional integrals on the right-hand side of (|334)) and ((335|) . 
Finally, since we are concerned on the leading term as iV — ^ oo in this article, assuming that we can deform 
the line segments in p.lSp or in p.29p to the real line, then we can once more interchange the order of 
integration and reconstruct the absolute value by integrating over r. 

In what follows, we say that the integral In as above satisfies the condition (vi): after an appropriate 
change of variable, the factor in Hi,{rS), depending on the variables w and the saddle points, becomes a 
real- valued variable (the line segments Li, coming from the integration in the neighborhood of the saddle 
points, should first be deformed to the real line). All the examples considered in the article satisfy this 
condition. The next corollaries immediately follow from Theorem 13.51 Remark 13.61 Theorem 13.71 and the 
use of the integral representation p.3ip . Similar results hold for the integral (|3.35p . 

Corollary 3.8. Under the foregoing assumptions (i)-(v) of Section \S.2\ and (vi) above, let 



N,n 



/ exp{-A^^p(t,)}|A(Or q{t)9{N^/^{t - to)) d"t 



where g(t) is analytic in C" and p{x) admits one saddle point xq of order /i — 1. Then, as N ^ oo. 



where 



^-nNp{xo) 
« n 

Ao= exp[-p('^Hxo)/^^lY,w^}g{w)\A{w)rd" 



Corollary 3.9. Under the foregoing assumptions (i)-(v) of Section \3.3\ and (vi) above, let 



lN.n = / cxp{-A^^p(t,)}|A(t)rg(t)d"i 



'(a.6)" ^.^1 

where p{x) admits two simple saddle points and Re{x'_|- — x^} > 0. Moreover, let p± = p"{x±) 

and Fy^TO be given in (|A.5|) . //Rc{p(x+)} — Re{p(a;_)}, then as N — > oo, 

(r, — T )'^™-^ p-mN{p(x+)+p{x-)) 
^N.2va l^„Jlii.,mj )r«+i^m(m- 1)/2 ^m+i/m(m-l)/2 



ASYMPTOTICS FOR PRODUCTS OF CHARACTERISTIC POLYNOMIALS 



17 



while 



1 — T yni{m-l) -mN(p(x+)+p(x-)) 



( ^ p^p_ y'n+iym{m-l)/2 ^(2m-l+iy(m-l)2)/2 

Wc conclude this section by applying the two last corollaries to the study of the asymptotic behavior 
of the generalized Airy function [TO] : 

^i'"^(^) = 7^ / o4"\^-,^t)dt (3.36) 

where it is assumed that the variables s = (si, . . . ,s„) E W\ Note that the generalized Airy function 
obviously reduces to the classical Airy function when n = 1. Moreover, for a = 1, the above n-dimensional 
integral is proportional to Kontsevich's matrix Airy function A(S), where S denotes a n x n hermitian 
matrix with eigenvalues si, . . . , s„ (see Section 4 in [35]). It is also worth noting that in the case where 
a ~ 1 and si = . . . = s„ = u, then the above integral representation can be reduced to a very simple 
determinantal formula (see for insante Section 4 in [TT] or Ref. 8 therein): 

Ai(2)(u,...,u) = (-l)"("-i)/2n!det ^ 



Ai(u) 



(3.37) 



Before displaying the asymptotic the generalized Airy function, we recall the following shorthand nota- 
tion: when A,B gM., {A + Bs) stands for {A + Bsi, ...A + Bs„). 

Proposition 3.10. Let x be a real positive variable. Then as x — > cio, 

Al (x + X s) ^ j-2^-jTi2(n+n(Ti-l)/a)/2 ^(n+n(n-l)/a)/4 ' (3.38) 

and for n = 2m 

Ai(")(-x + a;-i/2s) ~ (2™)i£2Z^^(2VJ)-™+"(™+i)/"e~*^'i(^\F^"^(TO/a;n/a;2is). (3.39) 

Proof. (|3.38[) and p.39p originate from integrals evaluated around one simple saddle point and two simple 
saddle points, respectively. For p.38p . set N ~ x^l"^ . Simple manipulations and the use of Proposition !^^ 
lead to 

/w-(n+n(n— 1)/q)/3 f 

Ai(")(7v2/3 + iV-i/3,s) = — / e^(^P^(*)/='+'P^(*»|A(t)|2/%J-(")(s;ii)dt. 

We thus have an integral like in CoroUarv 13.81 with p(<j) = —itj/3 — itj, q{t) = qJ^q"\s; it), and g{t) = 1. 
The function p has two simple saddle points at ±i. With xo ~ i, we have po = p"(xo)/2 = 1 which implies 
that the steepest descent path near xo would follow the horizonal line, as desired. We may thus apply 
Corollarv 13.81 to the case fi = 2, xq = i, Po = 1, p{xo) = 2/3, q{xo, ■ ■ . ,xo) = e"^^^*) and p.38p follows 
immediately. 

For dnSni), we also let N = x^/^. In the definition of Ai^"'(s), substitute tj by N^/Hj and apply 
Proposition [2?2l which yields 

)\r(Ti+Ti(ri— 1)/q)/3 r " , , 

(271-)" Ju" jJl 

where p{tj) = i (— tj/3 + tj). This function has 2 simple saddle points, namely x± = ±1. This time we 
have to consider both of them because they are already on the path of integration. We have p{x± ) = ±2i/3, 
p± = p"{x±) = =F2i. This means that the steepest descent path is given by 



^ = - 



|-l + re-"/4 : t€ (-oo, ^/2] | U |i + re^'^/^ : rG [-V2,oo)|. 
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By making the change of variables wj H> e^^^^^Vj (for saddle points x±) in p.29p . we see that the variables 
Vj follow the real line, so the assumption (vi) is fulfilled and p.39p follows from the previous corollary. □ 

4. Scaling limits 

In this section, we prove the theorems and corollaries given in the introduction. Most of the proofs rely 
on Corollaries l3.8l and 13.91 or similar results for integrals on the torus T", more precisely for (|3.35p . All the 
integrals considered here fulfill the assumptions (i) to (vi) given in the last section. Note that when we talk 
about deforming the contours of integration, we implicitly suppose that either the power v of the absolute 
value of Vandcrmonde determinant is even and the variables are real, or the integral representation p.31|) 
is being used. Also in this section is the assumption that a deformation of contours is made as long as the 
integrand is analytic and absolutely intcgrable over the concerned region of the complex plane. 

We will frequently use the symbol to denote a straight line path from a to b. Additionally, will 
denote a semi-circular path in the positive direction, starting at a and ending at b with radius \a — b\/2. 
Dyson index /? and its duality /?' := 4//3 will be used alternately. 

Before starting our computation we first review the integral representation of hypergeometric functions 
2F^°'\a,b; c] s) and iF^°'\a; c; s), due to Yan [51] and Forrester [21], especially [23l for more details. For 
a > 0, Re{:^i} > -1 and Re{i/2} > -1, 

2Fi^°^(a,6;c;5i,...,5„)= \ / , / i-F^"^ (a; t) W d"*, (4-1) 

where i^i = 6 — (n — 1)/q; — 1, 1^2 = c — b — [n — l)/a — 1 and 

n 

^..,...i/aW = n^r(i-i.r n i^'^-i.f^"- (4-2) 

i—1 l<i<j<n 

If Rc{!^i} > —1, the right-hand integral of (|4.ip can be analytically continued so that it is valid for 
Re{j^2} < ^1 but by replacing the interval [0, 1] with the counter-clockwise circle T, especially in the case 
of interest (a = -N) f4„21\: 



2Fi"\a,b;c;s) 



^i-K{h—c)r, 



Mn{b - c, c + 1 + (?i - l)/a, I /a) 



1 



1 J-^"^(a; s-l-t) Z?,,+(„_i)/„,,,a/a(t) (4.3) 



(2^i)" A" 

Note that the constant A/„(a, 6, a) is given in the appendix. Likewise, for Re{c — a} > —1, we have 

■ — -'17T71(1 1 /i 71 

.Ft\a;c+l + in-l)/a;s) ^ ^^-a.cM o.) X. 11 T - ^.)-|A(t)|^/%^^")(.; .) d". 

(4.4) 

4.1. Hermite /3-ensemble. It follows from the duality relation in Proposition 7 [10] that 



i^jv(si, . . . , s„) : = ( n n(^* - ■'*j))h«e 

= (-j)"^2^'"("-i)/4+n/2(p^,^)-igP.(.) / ^t'^^~'"\^(t)f^T'il''\-2is-t)<Pt. 



Set 

s 



s, 1-^ V2N(u + and t,- 1-^ V2iVt,-. 

pN 

By using Proposition 12.21 the weighted quantity 

(pjv(si, • ■ • , s„) := e-^P^^^'^KNisi, ...,s„) 
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becomes 

s 



^j^_i(V2N{u+—)) = (-zV2iV)"(^-'n2ViV)^"("-i)/2+"(r^,„)-ie"^" +^'=^"'/^''"^^/Ar,„, (4.5) 



where 



lN.n= / exp{-ivVp(t,)}|A(i)|^(i)d"t 

= cpJ\-^'-'(^J^jM-Nj2pmq{t)Hf,,{rAit))d-tyr, (4.6) 



and 



p(x-) =2i-2+4mx-lnx, g(t) = J| ij' o J'^^^'^'^ (-4is/p; t + m/2)- (4.7) 

Here I is a fixed integer, and it is assumed that I ~ for n = 2m while / = 0, 1 for n ~ 2m — 1. 

Since p'{x) = 4x + 4m — i, there are two simple saddle points in the bulk: x± = where 
u £ (—1, 1). At the rightmost soft edge (which corresponds to m = 1), these two points coincide and give 
a single saddle point xq = oi multiplicity 2. 

In the following lines, we proceed to compute the leading asymptotic terms in the bulk and at the soft 
edge. We start with the bulk of the spectrum. We set u = s'mO, 6 e (^f,^), so that x+ = ^e^'^, 
X- = ie*(^+'^) and 

p{x+) = -icos26'+ (1 + ln2) + i(6'+ isin26'), p{x_) = -icos26'+ (1 + ln2) - i{e + isin26' + 7r). 
It follows from 

P± ■^p"{x±) = 8e±"'cos6' 

that the angles of steepest descent are —6/2 at x+ and 6*72 at x_. Note that if a; = — |m + v, w G M, then 
the function 

1 13 

Re{p{x)} = 2^2 - - ln(w^ + -li^) + -u^ 

attains its minimum value at the point v = zt^\/l — . Therefore, as a possible path, we consider the 
straight line passing through —-^u and parallel to the real axis (together with irrelevant deformations at 
±oo). 

The segments of integration (—ooe~* 2^006^*2) and (— ooe*2 , ooe*2 ) can be deformed into the real axis 
near the saddle points. Since a;+ — x_ > 0, the assumption (vi) is fulfilled. Take p ~ — u^. According 
to Corollary EH for n = 2m (; = 0), 

lN,2rn (SA^)-'^'™^™-!)/^-™ g^p|_^^^2 _ ^^(^ + 2 In 2 - l7r)/2} (!^)/9'm(™+l)/2-m 

X C::){T,,„,fo4'^^\^ns■,r'\i-lr), (4.8) 
while for n = 2m — 1 , 

iNam^l - (8A)-^'("-l)'/2~„/2g^p|_^^^2 _ ^^^^ + 2 In 2 - in)/2} (Jllfim' ~l)/2^n/2 

X {'":;:')Tp,^^,Tp,^{-2^)(^-^-^^'[^^^^^ (4.9) 
where 



6n = Ar(26' + sin 261 + 7r)/2 + 6'(1 + (m - l)/3')/2, 9 = arcsinw. 

Hence, for n = 2m, 

^^{V2Niu+-^)) ^ ^'«,2,„7n^(/3') o-^,i'/^\»7rs;l("),(-l)(™)), (4.10) 

where 

*iv,2™ = (7rp)^'"("+i)/2-"7V^'™'/2exp{-mAr(l + ln2-lnA)} 
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and 

m J j-J-^r{l + /3'{m+j)/2) 



\ m J ■'-^ r(l + B'lm + 7)/2 



For n = 2m — 1, 



pN '"2iVcos0 

where 



X exp{-(2m- l)Ar(l + ln2 - ln7V)/2} ( VW2)"*^'""^^'(2i V'^p/2). 

We now turn attention to the soft edge of the spectrum (i.e., u = 1). In Eas. (|4.5p - (|4.7p . let u = 1, Z = 

and p = 2N-^/^. Then, 

p{x) = 2x^ + Aix - In X, q{t) = qF'^' ^\~2iN'^l^ s; t + i/2). 

At the double saddle point xq = — */2, we have 

p{xq) = 3/2 + In 2 + m/2, p"'{xo) = 16i. 

The angle of steepest descent are thus —5n/6 and — 7r/6. However, we see that Re{p{x + xq) — p{xo)} > 
for all X S M except at the origin. We thus choose ^ to follow the real line from — oo to oo. CoroUarv 13.81 
then implies 

lN,n ~ (8^)"'^'"("~l)/'^-"/3exp{-n^(3 + 21n2 + ^7r)/2}(27r)"Ai(2/^')(s). (4.13) 

Thus, 

^jv(V2iV(l + ^)) ^ i27ry\Tp,^r'Ai^'/^\s), (4.14) 

where 

= 7V'3'"("-i)/i2+"/6gxp{_^7V(l + In 2 - In TV - 2i7r)/2}. 

For (3 even, the scaling limits of the correlation functions for the H/3E immediately follow from (|4.10p 
and (|4.14l) . Let n = 2m = fc/3, we have the following relation 

Rk^Nixi,...,Xk) = Q^'^ n (^^J -2;;)^ [¥'Ar(si,...,Sn)]{,j^{^j. (4.15) 

One easily shows that as cxo, 

^2-T^yk/22mk+l)/i p-l3k/2^Y[l + ^/2))'=(2e)'^''^/27V-^'=^/2-'3fc(fc+l)/4-fc/2^ 



Gl3,N 



We consider the bulk scaling and let p = - \/l — u^. Some manipulations then lead to 



(^fRk,N{V2N{u+^)) ^ (/3/2)"^'=/2(r(l +/3/2))'=7„,(/3') |A(27rx) 1^0-^(1'/"'' (^^s; 1™, 
According to Gauss's multiplication formulas for the gamma function [5], 



nr(a+:^) =r''^+^(27r)^r(Za), ZeN, 



we have for I = /3/2, 



^ r/^v r(i + /3V2) .^/o^^ifeV2TT r(i + /3j72) 



iir(i + /3'(™ + j-)/2) ' ^A^i^r(i + /?(fc + j)/2) 
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This further iniphes 



i^f^^^^^^^"" + ;^)) - ^^^^^ \^{2^x)f,F^,'/^\^^s■, 1", {-ir){s}^{.} ■ (4.16) 
Note that the coefficient 

k—l 

:= (/3/2)^^-(^-i)/^(r(l + (3/2)r n ril + Pik+%2) ^^'^^^ 

is exactly same as that in the circular /3-ensenible (Proposition 13.2.3, 1121) as well as in the L/3E and the 
J/3E below. This strongly suggests the universality of &fe(/3). 
Similarly, for the soft edge, one sets p = 2N^^/^ and gets 



where 



i^fRk.N{V2Nil + ^)) ^ a,W) |A(x)|'^Ai(^/2)(5){,}^{.}, (4.18) 



-'=1-- --.^. v^v^ ■ 11 r(l + /3j72) 

4.2. Laguerre /3-ensemble. Based on the work of Kaneko [28|, it is easy to verify that if 

N n 

Kn{si, . . . , s„) = ( - Sj] 



LI3E 
1=1 j=l 



then 



K^is^, . . . , s„) = ^;+"'^-^ iff/^^ (-^; (2//3)(Ai + n); si, . . . , s„) , (4.20) 

for instance, see Proposition 13.2.5, [23]. By making the change Si i— >■ s^/A^ and taking the limit N — ?> oo, 
one readily proves the first formula of Theorem 1 1.31 Note that this result could be obtained by taking the 
integral duality formula in [TO] and following the asymptotic method developed in the previous section for 
the case where p{x) admits a simple saddle point. 
By using (|4.4p . we get the following integral formula 

« n 

K^{s^, s„) = ^-" J^^ n - t,r-'+^''^^^+'^/'\A{t)fo4'^'"> {s; t) d^t (4.21) 



where 



(27r)-"e''^"^ Wx 



A/„(Ar,-l + /3'(Ai + l)/2,/3'/2) Wx,.p,N ' ^' ' 



Now let 

l<j<n 

The application of Proposition 12.21 then yields 



where 



^^_,(47V(^+^))=AA._,(4A^)"'^/^e-2"^" n {^+^)''^'lN,n, (4.23) 



/A..n= / CXp{- Ny^p{t,)}\ Ait) fqit)rt 

= ^"^'"M / cxp{-A^^p(t,)}g(t)iJ^,(rA(t)n(^i.)"*""'^/')d"tWr, (4.24) 
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and 

n 

p{x) =\nx- ln(l -x)- Aux, q{t) = ^-^(l - i^.)-'-i+'9'(Ai+i)/2 ^j^i^f^'^^s/p; t - 1/2). (4.25) 

3 = 1 

Since p' (x) = ~ 4ii, there are two simple saddle points in the bulk, namely x± = ^ zt 

with u € (0, 1). By letting u — 1, which corresponds to the soft edge, we find that the two saddle points 
become one double saddle point xq = 

We first focus on the bulk of the spectrum. Set u = cos^ 9, 9 <E (0, f ), hence x± = „ "'" ^ e^'^, and 



p{x±) = -2cos^9±i{29-sm29), p± p"{x±) = ±16iu^,J^. 

We see that the angles of steepest descent are |7r at x+ and jn at x-, so we choose the following path of 
integration: 

«,l/(2cose)| I ^-l/(2cose)| I ^l/(2cose) I I ^1 
■^1 U l/(2cose) U -l/(2cose)U l/(2cose)- 

Actually, when x G jf^V(2'=°'^*) or -Sfjy^jcosS)' have 

Re{p(a;) - p(x+)} = In r - Aux + 2u > 0. 

\l~x\ 

On the semicircle . we can write x = „ ^ „ e"''' with € [0,7rl. Since 

l/(2cos(?) ' 2cost^ ^ L 7 J 

(/(cos 0) := Rc{p(a;) — p(x+)} = — ^ ln(l + 4 cos^ 6' — 4 cos 9 cos 0) + 2(cos 9 — cos 0) cos 6*, 
it follows from 

,, 8(cosfl!) — COS0) cos^ 6* 

q (cos (7)) = 7, 

l + 4cos2 6'-4cos6icos0 

that g{cos(f)) attains its minimum value at cos0 = cos0. Similar results hold for Rc{p(.t) — p(x^)}. 
Case 1: n = 2m. 

The line segments of integration (— cxje'i'^, ooe'*'^) and (— cx^e' ooe' *^) become the real line by making 
the following change of variables: Wj i— >■ WjC^i'" and Wj >• Wje^^'^ , respectively near x^ and We take 

p = fy^^^ and / = 0. Then, according to Theorem 13.71 we have 

/jY 2 _/V~'3'™(™~l)/2-™g4mAfM^2y^^^/3'™-^i ^ '^P ^ff 'Tra(m+1) /2-m 

/2m\/-p n2 x-(2//3') 



X (^™)(r^'.™)%-^r^'(*'rs;l",(~l)™). (4.26) 



; From 



we also get 



Hence, 



where 



Wx,+n.AN ^ f ,^.nN TT ^ + ^1 + n + P.] /2) ^ ^ TjN + P\X^ + j) /2) 
Wx,,f>,N ^"^ /i r(l + Ai+/3j72) -y r(/?'(Ai+j)/2) ^ 



4 (O^^-nr.nN fV r(l + /372)r(l + TV + /3^(n - j)/2) 

^"^^'^^ ^ n wnm 

- (r^.„)-i7V'''"("-i)/'*+"/2 cxp{-niV(l - In AT _ ^tt)}. (4.27) 



*iv,2„. ^ (7rp/2)'^''"("+l)/2-m^/3'm(m+Ai)/2g^p|_2„iiV(l _ IniV)}. 
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Case 2: n = 2m — 1. According to Theorem 13. 7[ two types of integration domains lead to a significant 
contribution as iV — > oo. They correspond to the intervals of integration 

(— ooe ,ooe ) x [—ooe * , ooe ■> j , (— ooe * ,ooe " j x (— ooe * ,ooe j 

which can be deformed to (whenever m > 1) 

(-ooe*3",ooe*3")"-i x (-ooei"+^(^-2»\ ooe^'+^^^^^e))™ 

and 

Then, by changing the variables as in n = 2m case, we get 
where 

On = iV(26l-sin26l) + /3'(Ai + 1)6/2 + I3'{m - l)(6i - 7r/4) + 7r/4, 6* = arccos^u- 
One finally obtains 

^^._,(4iV(.+ ^)) ^ (e'(^""^'^)i;i^/fH7r.) - (z..)), (4.30) 

where 

= ('r')r/j',m~ir0',™(r/3',2™-i)-^(7rp/2)^'('"'-i)/2-"+iV2(2V^)-^'/2+i 

^ ^/3'm(m-l)/2+/3'(2m-l)Ai/4 exp{-(2m - 1)7V(1 - hi TV - ITt)} ( _iV) " f^™- 1)' . 

We now consider the soft edge of the spectrum (u = 1). In Eqs. (|4.23p - (|4.25p . let w = I, ^ = and 
p = 2(2iV)-i/3, This yields 

n 

p{x) =\nx- ln(I - x) - Ax, q{t) = z"" [| tj^(l ~ i^.)-i+/^'(Ai+i)/2 oJ-^2//3')(^^^^. ^ _ ^^g). 

At the double saddle point xq — 1/2, we have 

p{xo) = -2, p"\xo) = 32. 
Thus, the angles of steepest descent are 27r/3 and 47r/3. Since on the circle x = ^e**^, we have 

Rc{p{x + xo) — p(xo)} = — In ^5 — 4 cos — 2 cos + 2 > 

whenever </> g (0, 27r). We choose to be the following path: it starts at 1, arrives at l/2+i0+ by following 
a straight line, along the centered circle of radius 1/2 counterclockwisely, then follows a straight line from 
1/2 — iO+ to 1. Applying Theorem 13.51 and Remark 13.61 noting that the line segments of integration near 
the saddle point can be chosen to be (— cxji, ooi), after the change of variables: Wj l%~^l^iwj one obtains 

/JV,„ - 2-^'"("-l)/S-'^'"(^i + l)/2+2n/3^-/3'n(n-l)/6-n/3 cxp{2nAf} {2-kY ^X^"^ ' '^'\s) . (4.31) 

Thus, 

(^^(4Ar(i + j-)) ^ (2^)"(r^,„)-'Ai(2/^''(s), (4.32) 

where 

= 2-^'"("-l)/^-'''"/2+2"/3^^'n(n-l)/12+/3'«Ai/4+n/6g^p|_j^^(^ - In - m)} . 
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When /3 is even, the scaling hmit of the correlation functions for the L/?E immediately follow from (|4.28|) 



and (|4.32l) . For the bulk case, let n = 2m = kjS and p — niaking use of 

RkM[xi,...,Xk) = TT^^''^^^ n 'y^o-^if['PN[si,...,Sn)\,,^,y (4.33) 

and 

^^i-/^'^ _ {2ti)-^ {13/ 2)- ^^/'^{T{1 + _0/2))fce'''=^7V-'^''^-'"='/2-(^i+i)'', 
one can show that 

(^)'i?fc,^(4iV(u + ^)) ^ \A{2nx)f^Fi^'^\ins- 1™, (4.34) 

For the soft edge, one lets p = 2{2N)^^/^ and gets 

(^)'i?,,^(47V(l + -^)) ~ a,(/3) |A(x)|^Ai(^/2)(s){,}^{,}. (4.35) 
Notice that the coefficients afe(/?) and &fc(/3) are the same as those in the H/?E. 

4.3. Jacobi /3-ensemble. The J/3E case is very similar to the L/3E. First, Kaneko [22 proved that 

N n 

- ^"iV"/';/o('^ (-^' + A, + n + 1) + ^ - 1; (2//^)(A, + n)- s) . (4.36) 

By making the change Si i— > Si/N'^ and taking the limit iV — >■ cxj, one readily proves the second formula of 
Theorem O 

By using (|4.3p we get the following integral formula 



i^^(s) = i3wz-" / ni-^'(^^+i)/2-~(l-t,0^-'+^'(^^+^^+')/'|A(i)|^\j-('/^')(-iV;s;l-t)d^^ (4.37) 
where 

^ (27r)-"e-(^'"(^^+i)/2+»(^-i)) 5A.(Ai+n,A2,/3/2) 

7\/„(/3'(A2 + l)/2 + 7V-l,-l + /3'(Ai + l)/2,/372) 5jv(Ai, A2, /3/2) ' ^ 
For the weighted quantity 

l<j<n 

application of Proposition 12.21 gives 

^N-^ i-+^)= Bn-^ \{ + ^t'" (1 - " - ^t'" (4-39) 



where 

/iv^n = cp, T-^'-' y exp{-7V}_Jp(i,)} q{t) Hp,{rA{t) \ \ {it,)-^"-'>r^) d"t ]dr, (4.40) 

and p{x) = Inx — ln(l — x) — ln(l — u -\- ux), 



/ exp{-7V^p(i,)}g(t)ff^,(rA(i)n(^*.)"*""'^^')'^"0^^' 



1 + 

qit) = IT i-MA2 + l)/2+Z(^ _ , y(A,+A.+2)/2-;-2 ^_^(2//?') ^. ^ ^ ). (4.4I) 

piV 1 — u + ut 



Since 
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there are two simple saddle points x+ = ^/(l — u)/ue™^''^ and X- — 1/(1 — u)juc''^'"l'^ in the bulk with 
we (0,1). 

Set u = cos^ 6*, 6* € (0, 7r/2), then 



This allows us to take the angles of steepest descent jtt — 9 at x+ and 6* — -jTr at x_. We choose the path 
of integration: 

■^1 U -^t^nr' U -^-Tle U ^ta, 



Actually, when a; € Jf^*^^"^ or ^tinC" '^^ see that 



^rn cr,n + 

tanS' 



Re{p(x) — p(a;±)} = In ■; — ln(l — u + ux) > 0, 

\l-x\ 

while on the circle {x : \x\ = tan6'}, setting x = ta.nd e'-'^ with (f> G [0, 27r), 

Re{p{x) - pix±)} = ln(l - 4m(1 - u) cos^ 0) 

attains its minimum at = ±7r/2. 

Let p = — . , notice the polynomial 

^u(l-u) 

^ pN 1 — u + ut' ^ p 1 — u + ut' 

and 

Bn - (r^, „)-l2-^'"(^i+^^+2)/2-^'"("-l)/4-"(2W-3/2)^/J'»(n-l)/4+n/2 exp{j^(/3'n(A2 + l)/2 + n{N - 1))}, 

we are ready to compute the bulk scaling. 
For n — 2m, 

Vn{u +^)- *«.-7m(/3') o4'^^\^^s-^ 1". (-1)'"). (4-42) 

where 

^N2m. = (7rp)/3'™(™+l)/2-™_/V'^'™V22-/3'™V2-/3'm(Ai+A2 + l)+2m(l-2iV) 

while for n = 2m — 1 

^-^-'("+ ^) - <,„^(e'^"+^(^-^^'i?i^/f .) - e-««-(f-^)'4?^'',\z,r.)), (4.43) 

where 

C!.™-. = ('r')^/3^™-l^/3',™(^;3',2™-l)'^(^p)^'("-^)'/'+(^'-''"+^)/'y2(2^AI)-^'/2+l 

^ ^/3'm(m-l)/2j-l^_-[^^n(Af-l)2-/3'm(m+l)/2-/3'(2m-^l)(Ai+A2)/2+(2m-l)(l-2W+20 + l(|2^ _ w)"'/^ 

Here 

6)^, = 2{N~1)9+{1 + P'{m + X2)){e ~ Tr/4) + l3'{Xi - A2)6'/2, 9 = arccosv^. 

As previously, when j3 is even, scaling limits of correlation functions in the bulk for the J/3E immediately 
follow from ()4.42p . Let n = 2m = fc/3, we have 

(fc + A^)! 5jv(Ai,A2,/3/2) -pr . 
i?,,^(xi, . . . , xfc) = AT! ^,+^(Ai,A2,/3/2) ^ JJ^,^"^^' " ""'^ f^^^'^' ' " ' ' '"^^ {^}-^->- ^ ^ ^ 

Notice 

■S'Af(Ai, A2,^/2) ^ 7r-fc(p(i _^ ^^2))^2''(^^~^+'')''+^(^i+^"+^)''(^A^)~''''/^, 
5'fc+Ar(Ai, A2, ^/2) 

in the bulk taking p = — , ^ , one obtains 

^ Vm(I-m) 

{^t^KN {^+-^)- b,il3) \A{27:x)\^4'/'\^7:s■, 1", (-1)"){,}^{.}. (4.45) 



26 



PATRICK DESROSIERS AND DANG-ZHENG LIU 



The use of Corollary 12.31 finally allows us to rewrite the right-hand side in terms of a iFi hypergeometric 
function as in Corollarv ll.51 

We haye thus completed the proof of Theorems 11.11 II. 2[ 11.31 and Corollaries 11.41 11.51 stated in the 
introduction. 

4.4. The universal multivariate "kernel" in the bulk. Let us conclude this section by exhibiting a 
uniyersal pattern observed in the bulk of the classical ^-ensembles when n = 2m — 1 is odd. From the 
asymptotic results contained in Eqs (|4.12p . (|4.30p and ()4.43| ). one easily establishes the following theorem. 

Theorem 4.1 ("Kernel" in the bulk). Assume that n = 2m — 1 is odd. Let A = V^N, 4:N, I and 
p = — u^, iy^^^ "'^'^ ^ ^ — J for the H/3E, L(3E, J/3E, respectively. Moreover, let s — (si, . . . , s„) 
and t = (ti, . . . Then as N — > oo, 

- ^ (-^) ( — t) - eTI (-0 (— ^) } (4.46) 

where E^'^ denotes the generalized exponential defined in (j2.19p and ^'^ 2,71-1 given in (|A.9[) . 

The last theorem obviously generalizes the standard result valid for the unitary (i.e., /3 = 2) ensembles 
and according to which the polynomial kernefl K]\[{x,y) asymptotically tends to the sine kernel 
when it is rescaled in the bulk of the spectrum (see for instance [23ll41) ). 

5. PDES AT THE EDGES AND IN THE BULK 

Since Kaneko and Yan [HI], we know that the hypergeometric function qF^^^^\{2/ l3){Xi + n);s) 
satisfies the following holonomic system of PDEs: 

d^F 2.^ ^ .dF ^ 2 ^ I f OF dF\ „ , , 



Now, we also know from Theorem 11.31 and Eq. (|1.7p that the expectation value (^at, when rescaled at the 
hard edge of the L/3E or J/3E, is given by (si ■ ■ ■ Sn)^^^ ^ qF[^ ^'^\{2/ + n); — s). Consequently, the 
limit of the rescaled ip^ satisfies the following system of PDEs: 

d^F 2dF Ai,Ai+2 2^ 1 / dF dF\ „ , , 



Sfe 



dsl /3 dsk \ 13 13 skj P sk - Sj \ dsk dsj 



Similar results hold in the bulk and at the soft-edge of the classical /3-ensembles. This was first shown 
in by exploiting Kaneko's system of PDEs for the hypergeometric function 2^1 [IH] 0. More explicitly, 
by properly rescaling tpN and then taking the limit — > 00, one arrives at the conclusion that the limit 
of the rescaled function ip^ satisfies in the bulk (n = 2m), 

d^F ^ 2 ^ I fdF dF\ „ , , 

while at the edge, it satisfies 

d'^F ^ 2 ^ 1 f dF dF 



dsi "'^^ ' /? . ^ Sk-Sj \dsk dsj 



IpoIIowing our notation, Kn{x, y) is equal to W-^"(;/)^iv-i(^) ^.^j^ „ ^ ^ 

"^It was assumed in 1361 that th 
make this assumption superfluous 



"^It was assumed in 1361 that the scaling limits of (/jjy should exist in the bulk and at the edge. Theorems ll.2l and ll.ll now 
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By summing up n equations in (|5.4p . one easily shows that the multivariate Airy function is one possible 
solution of the following (single) PDE: 

DoF = pi{s)F. (5.5) 

This PDE was first given [TUl Eq. (5.17)] as an equation satisfied by the Airy function defined in (jl.l2p . 

In the one-dimensional case (j5.3p and (j5.4p respectively reduce to the well-known differential equations 
satisfied by the sine (or cosine) and Airy functions. However, in the higher-dimensional case, it seems 
difficult to find all solutions of (|5.3p and (|5.4|) . one of which being our limiting expectation value. 

Let us consider the limiting bulk expectation value with n = 2. In this case, we can express F{si, S2) = 
oJ-Q^^^^(i, —i; si, S2) in a more explicit form. Firstly, according to (|2.12|) . we have 

^^(si,S2) = e-*^i+'^^oJ"f/"(2i,0;si-S2,0), 

which only depends on si — S2. Set F{si,S2) — /(si — S2), so that f{x) is an analytic function with 
/(O) ~ 1 and f{x) — f{~x). Secondly, (|5.3p implies that f{x) satisfies 

/" + M)/' + /^0, (5.6) 

X 

which can be reduced to the Besscl equation (cf. (4.5.9) in [2]). From this, we get 

/(x)=2l-5r(| + i)x^-|j|_i(x), 

where Ja{x) is the Bcssel function of first kind of order a. Furthermore, 

oJ-^^/'^(*, si, S2) = 2f -^r(| -f i) (si - S2)^-^ J|_ 1 (si - S2). (5.7) 

The right-hand side of (|5.7p has been first obtained by Aomoto [3] (at zero) and by Su [47] in the bulk of 
the H/3E with < /3 < 4. 

The case n = 2 at the soft edge has also been previously studied. Indeed, Su [l^ has proved for the 
H/3E, with the aid of Dumitriu and Edelman's tri-diagonal matrix model, that the scaling limit of lpn{s) 
has a single integral representation: 

1 rl + joO gTi2-^-^(si+S2)-t7(«l-'S2) 



4,^3/2^ 



2,1 dz. 

2/3 + 2 



This integral was first defined by Kosters |35j ; it is believed to be proportional to our 2-dimensional integral 
representation. 

We end this section with a few remarks on the (3 = 00 case. It is well known that the parameter /3 of 
Random Matrix Theory can be interpreted in Statistical Mechanics as the inverse temperature of a log-gas 
system. Thus /3 = 00 corresponds to the completely frozen state of the system, which is the state where 
the particles no longer move. Moreover, given that Pi°°-'(s) = 771^(3), we have 

= ? m.(l») • ^'-'^ 

We claim that 

o-F^)(y;.) = ^ ^ n^™'- (5-9) 

■ <Tes„i=i 

Actually, the latter identity follows from the expansion of the right-hand side of (|5.9p , the use of (|5.8p and 
the following easily established fact 



n 
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By using ()5.9p . we easily get 

n 

Ai(-)(s) = nAi(«^)' (5-10) 

where Ai(x) denotes the one- variable Airy function of first kind (Bi(.T) for the second kind). This result is of 
course consistent with the physical phenomenon of complete separation of the particles at zero temperature. 
However, the situation is a little more complicated in the bulk. For instance, when n = 2, the use of (|5.7p 
implies that 



oJ"o°°''(i,-i;si,S2) = 2 2r(^)Vsi~^ - S2) = cos(.si - 52), (5.11) 

in which the variables cannot separate. We refer to jl4j for more information on large /? asymptotics. 

Another way of putting this is that when /3 = 00, the systems of PDEs such as (|5.3|) and (|5.4p can be 
directly solved. Indeed, the linearly independent symmetric solutions are respectively 

1 y e'("-(i'+-+"-«'-"-u+i) — =o-F^°°^(l^(-l)"-^';^s), .7 = 0, 1, . . . , n, (5.12) 



,1, 



and 



-7 X! ^Ks<T(i)) • • • Ai(s^(j))Bi(s^(j+i)) • • • Bi(s^(„)), j 0, 1, . . . , n. (5.13) 



Note that the equality in (j5.12p comes from (|5.9p . 

Now a natural question arises: Can we guess the similar results to (j5.12p and (j5.13p for general /3? In 
particular, are o-T^o'^^^' ((— l)"* , l""-* ; is) where j = 0, 1, . . . , n all linearly independent symmetric solutions 
of JO])? Note that if j = or n then 

oJ-f /^'((-l)^ l"-^";i.s) ^ e±*fi('*) 

satisfies (|5.3p . Moreover, according to our bulk scaling at zero for the H/3E and for n = 2m or 2m — 1 

oJ-f /2'((-l)™,l"-™;±is) 

are also solutions of (|5.3p . see (|4.10p and (|4.12p (0 = 0, Z = 0, 1). So we give the following 



Conjecture: The set of symmetric solutions of the system of PDEs (|5.3p is spanned by n + 1 linearly 
independent functions 

oJ■f/''((-l)^l"-^■;^s), .? = 0,l,...,n. 
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Appendix A. Notation and constants 
Most of the constants used in the article can be derived from Selberg's integral [23ll4T| : 

N 

2X3 



5iv(Ai,A2,A3) : - / f\x^'{l~x,)^^ FT \x, ~ x,\ 

•^[0.1]" ^=l 1<^<J<N 



^ W r(i + A3 + jA3)r(i + Ai + .7A3)r(i + A2 + jh) 

7=0 + ^3)r(2 + Ai + A2 + (A^ + j - 1)A3) ■ ^ ' ' 

In articular, one readily shows that 
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and 

G,,A. = /?-^/-/3^(^-i)/4(2.)A^/2 g ^(^+|^/^' + ^|^/') . (A.3) 

For the Gaussian case, it is often more convenient to use the following integral: 
/■ " 1 

/ n^"'"'^' n l^.-x,f rf^i---rfx„= ^„^^„^„_^)^,^^^ r,,„, Re{z}>0, (A.4) 

where 

The Morris normalization constant is 

^ r(l + a)r(l + a + ja)r(l + 5 + ja) ^ ' 

The constants for the the soft edge are 

_ fiV'3'"("-i)/i2W6gxp{-nAf(l + ln2-lniV-2i7r)/2} H/3E 

~ I 2-/3'"("-l)/6-^i'«/2+2«/3^/3'n(n-l)/12+/3'nAi/4+n/6g^p|_^^^ - In iV - ITt)} L/3E 

In the bulk with n = 2m, we have 

{(•^p)/3'm(m+l)/2-m^/3W/2gj^p|_^jY(l + ln2 - IniV)} H/3E 
(7rp/2)'9™('"+i)/2-"W3'™("+^i)/2exp{-2mA^(l - IniV)} L/?E (A.8) 
^^p-^P'm{m+l)/2~m]\jl3'm^ /2-l3'm{Xi+X2 + l)+2m{l-2N) j^g 

where p = — w^, ^\/^^, h , ^ , respectively correspond to the H/3E, L/3E and J/3E. 
The constants for the bulk with n ~ 2m — 1 are 

f ('r')r^>,_ir^,™(r^,2„._i)-i(^p)'''(™'-i)/2-(2'»-i)/2iv'5''»(™"i)/2 

xexp{-(2m- l)Ar(l + ln2 - In Ar)/2} (y]v72)-(2™-i)'(2i^7rp/2) H/?E 

('r')^/3',m-l^/J',™(^^',2™-l)-l(W2)^'(™^-l)/2-™+l^/2(2V^)"/^'/2+l 

x]v'3'™("^-i)/2+/3'(2m-i)Ai/4gxp|_(-2^„ _ IniV- i7r)} (_7V)-(2™-i)i L/3E 

(''r')r/?^™-ir^i^™(r/.^2™-i)-H^p)^''™-'^'/'+'^'-'™+'^/'V2(2V^)-/^V2+i^^iv„(™ 

Xj-l(_l)n(JV-l)2-/3'm(m+l)/2-^'(2m-l)(Ai+A2)/2+(2m-l)(l-2Af+2/) + l(--[^ _ i^ J/3E 

(A.9) 

For the hard edge, the constants are 



_ r ^^^±^ L/3E 

^^'.^ = ) S»(A;+n,A2;/3/2) (^'lO) 
I Sjv(Ai,A2;/3/2) -^P-^' 



Finally, the universal coefHcients are 

k — 1 

bm = (/3/2)^'^('=-i)/^(r(i + p/2)r n r(iT/3V+ J j/2) ' ^^-^^^ 



and 
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